Quantum phases in a doped Mott insulator on the Shastry-Sutherland lattice 
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We propose the projected BCS wave function as the ground state for the doped Mott insulator 
SrCu2(BOg)2 on the Shastry-Sutherland lattice. At half filling this wave function yields the exact 
ground state. Adding mobile charge carriers, we find a strong asymmetry between electron and hole 
doping. Upon electron doping an unusual metal with strong valence bond correlations forms. Hole 
doped systems are d-wave RVB superconductors in which superconductivity is strongly enhanced 
by the emergence of inhomogeneous plaquette boi 



The physics of doped Mott insulators is at the heart 
of some of the most interesting and challenging problems 
in condensed matter physics In the extreme "atomic" 
limit, the localization of electrons in the Mott insulat- 
ing state leads to a macroscopic degeneracy of their spin 
or orbital degrees of freedom. Lifting this degeneracy 
via non-local quantum fluctuations gives rise to the rich 
variety of competing ground states in real materials. In 
many cases magnetic or orbital order prevails. In systems 
with geometric frustration and strong quantum fluctua- 
tions the possibility of a ground state without such order 
exists. Strong quantum mechanical fluctuations are also 
believed to induce unconventional superconductivity in 
doped Mott insulators and are related to the emergence 
of resonating valence bond (RVB) type fluctuations 0] ■ 
The investigation of doped Mott insulators without long 
range magnetic order at half filling is therefore the most 
direct way to uncover the nature of the strong local quan- 
tum fluctuations and their impact on the charge carrier 
dynamics, including superconductivity. It allows one to 
address questions such as: Are all doped RVB states 
superconductors? What is the role of spatial inhomo- 
geneities? or What determines the stability of states such 
as the d-wave resonating valence bond?. By analyzing a 
specific example, doped SrCu2 (603)2 with Cu-spins on 
a Shastry-Sutherland lattice, we give answers to these 
interesting questions. 

SrCu2 (603)2 is a Mott insulator with a low tempera- 
ture spin gap As ~ 30K[1I3]. The Cu'^+ states with s = 

spins interact predominantly within two dimensional 
CuBOs-planes and form dimers along the bond with 
dominant spin-spin interaction. The topology of the Cu- 
sublattice is shown in Fig.[T] The dominant, intra-dimer 
exchange interaction J' is along the diagonal bonds, while 
the smaller inter-dimer exchange J is along the x- and 
y-axes. There are two remarkable aspects of this ma- 
terial. First, the exact ground state of the Heisenberg 
model, on the lattice shown in Fig. [T] was obtained by 
Shastry and Sutherland^]. For sufficiently large J' the 
exact ground state energy per site is Eq/N — — | J'. The 
ground state wave function is an orthogonal dimer state 
|$od), i.e. a direct product of singlets along the bonds 
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correlated metal correlated superconductor 
with localized singlets with plaquette bond order 

FIG. 1: Shastry Sutherland lattice with efi^ective hopping Xii 
(for electron doping) and local pairing strength Aij (for hole 
doping). The width of each line determines the magnitude 
Xij or Aij , respectively. For electron doping, strong correla- 
tions along the diagonals occur, reminiscent of the undoped 
valence bond crystal. For hole doping, the nature of the pair- 
ing state changes, leading to d-wave superconductivity and an 
inhomogeneous plaquette pattern of the bond strength. Pla- 
quette states with large Aij around the two different diagonal 
orientations are degenerate. 



connected by J': |<i>od) = life) (iTi)., - UT).,) /\/2. 
Shastry and Sutherland demonstrated that this is the 
ground state for J' j J > 2. More recent numerical 
results d, 0, H, 0, [13, m support a transition for J' /J 
between 1.42 and 1.47. In the state |3>od) singlets are lo- 
calized, making this Mott insulator a valence bond crys- 
tal without long range magnetic order. The second in- 
teresting aspect of SrCu2(B03)2 is that J'/J ~ 1.57 
{,]' ~ 85K and J ~ 54K|4|]) is close to the regime where 
|$od) ceases to be the ground state. The valence bond 
crystal is comparatively fragile and competing states 
become relevant. Interesting proposals for competing 
phases at half filling include helical ordered states IJ, 3 , 
plaquette singlet states[l,|Tl|,[l3| and a deconfined spinon 
phase [ist. 

In this paper we investigate the stability of the va- 
lence bond crystal, and the emergence of superconduc- 
tivity, flux phases and inhomogeneous states in doped 



SrCu2 (603)2- We use a resonating valence bond vari- 
ational wave function that includes the exact ground 
state energy for the half filled parent compound, mak- 
ing the approach particularly well controlled for small 
doping. We demonstrate that upon adding charge car- 
riers, the system displays a strong asymmetry between 
electron and hole doping. For electron doping a non- 
superconducting but metallic phase with strong valence 
bond correlations forms, while an unconventional super- 
conductor emerges for hole doping. In the latter case 
a spontaneous breaking of the spatial symmetry occurs, 
where the paring gap forms a plaquette pattern, while 
the charge density remains uniform (see Fig.[T]). We then 
show that for hole doping superconductivity of this in- 
homogeneous state is strongly enhanced compared to the 
homogeneous state. 

In order to describe doped SrCu2(B03)2 we start from 
the t-J model [ll: 



H 



(1) 



where X^o-cLciCT < 1, the tilde indicates that these op- 
erate in the projected Hilbert space, and the remaining 
notation is standard. In addition to EqUwe also include 
three site hopping terms (not shown explicitly) that 
are of same order in t/U and become relevant for larger 
doping. We include hopping elements t and t' for the 
bonds with exchange interaction J and J', respectively. 

Hole doping in the t-J model is achieved by taking 
out electrons. To describe electron doping we perform 
a particle-hole transformation and take out holes. This 
transformation changes the sign of t' but not oft. Thus, 
one sign of t' corresponds to hole, the other sign to elec- 
tron doping. Since t' is the largest hopping element we 
expect an asymmetry between electron and hole doping. 
The effect is already visible if one analyses the single hole 
state on a 4-site lattice. A doped carrier is localized on 
the diagonals for t' < and delocalizes for t' > 0. Exper- 
iments so far determine J'/ J = [f /t)^ but not the sign 
of t'. 

In order to determine the crucial sign of the diagonal 
hopping elements t' we perform electronic structure cal- 
culations for SrCu2 (603)2- cuprate superconductors 
such calculations give the proper Fermi surface and band 
dispersion despite their problems with the Mott insu- 
lating state and mass renormalization effects. We assume 
the same is true in our case. In Fig. [2] we show our results 
in comparison with a tight binding model where we ig- 
nore correlation effects, making the system a (semi)metal 
instead of a Mott insulator. In the LDA calculation we 
ignore magnetic order. Except for the detailed behav- 
ior around the M-point the tight binding model gives a 
very reasonable description of the electronic structure. 
It also demonstrates that other bands, not included in 
the tight binding approach, are well separated in en- 
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FIG. 2: Bandstructure of SrCu2 (603)2 obtained using den- 
sity functional theory in comparison with a tight binding fit. 
Only the in-plane dispersion of the bands is shown. The bands 
close to -Ef are made of Cu-3ci^2_y2 and 0-2p states, while 
the well separated bands around — 1.25eV consist of 0-2p 
and Cu-3da;z and Sd^z states. 



ergy. A fit close to the F-point yields t = 0.09eV and 
t' = 0.104eV > 0. The ratio t' jt ~ 1.15 is slightly 
smaller than the value 1.25, obtained from the ratio of 
the exchange constants. However, a fit where one forces 
t' jt = 1.25 still gives a very good description of the LDA 
band-structure. In what follows we measure energies in 
units of t, use t' — ±1.25t for hole and electron doping 
and Jij = 0.3\tij\. The qualitative behavior of our re- 
sults is unchanged if we use the experimental values for 
Jij and the unrenormalized values for tij listed above. 

We now include the effects of strong correlation into 
the calculation by using a resona ting valence bond 
ground state wave function [Tsl 17, 3| 



I*) =P|$BCS), 



(2) 



where P projects out all doubly occupied states. |$bcs) 
is the BCS-ground state wave function. The RVB de- 
scription is particularly appealing as it reproduces the 
exact ground state at half filling (see below). I'I'bcs) is 
the fixed particle projection of the ground state of 
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The effective hopping elements Xij and pairing gaps Aj^ 
are the variational parameters of our many body wave 
function and are determined by minimizing the energy 
£; = (^f \H\ E is evaluated using the Monte 

Carlo approach of Refflil for 12 x 12 up to 18 x 18 lat- 
tices. The many body wave function Eql2]also allows de- 
termination of the long range correlations between pairs, 
enabling us to distinguish between local pairing and off 
diagonal long range order (ODLRO). Only the latter cor- 
responds to superconductivity, yields a Meissner effect 
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FIG. 3; Superconducting order parameter, signaling ODLR 
for hole doped systems. Even though the wave function for 
electron doped systems is characterizd by a pairing strength 
7^ in the BCS wave function, no superconducting order 
exists. From the analysis of the low energy Drude weight we 
conclude that the system is a metal. For hole doped systems 
d-wave superconductivity emerges. The superconducting or- 
der parameter is strongly enhanced in the inhomogeneous pla- 
quette phase. The inset shows the local spin correlations for 
electron and hole doping. The thickness of the line refers to 
the magnitude of (SoSi). All nearest and next nearest neigh- 
bor correlations are antiferromagnetic. 

and flux quantization. To this end we analyze the pair- 
ing correlation function, i^a.b (r ^ r') — |&J.a^r',b| 
with 6r,a — \ (cr+a,j,Cr| — Cr-i-a,TCr.i) and a, b nearest 
neighbor vectors. |'0s|^ « i^a,b (R = r — r' ^ oo) deter- 
mines the superconducting order parameter 

Undoped system: At half filling, the variational wave 
function in Eq[2] yields the exactly known ground state 
energy Eq with an accuracy 10""^. We further obtain 
the expected spin correlations (S^ • Sj) which equals — | 
for sites i and j on the same dimer bond and vanishes 
otherwise with error bars of lO""*. Thus, at half filling the 
RVB wave function, Eq[21 reproduces the exact ground 
state of the Shastry-Sutherland model, making it an ideal 
starting point to investigate doped systems. 

Doped valence bond crystal: At finite electron doping, 
we find that the renormalized diagonal hopping x' that 
enters |$bcs) as a variational parameter is significantly 
enhanced compared to the bare values. For doping values 
X ~ 0.1 we find x'/x — l-5i'/^i see also the left panel of 
Fig. [TJ The spin correlations (S; • Sj) shown in the inset 
of Fig. [31 are large for sites i and j on the same diagonal 
bond and very small otherwise and mimic the behavior 
of the orthorgonal dimer state of the undoped parent 
compound. 

The pairing gap Ay is finite (for example, at electron 
doping a; = 0.1, we obtain the gaps along diagonal bonds 
^i,i+x±y ~ O.St ). It is remarkable that this many body 
state with finite local pairing gap does not show any indi- 
cation for off diagonal long range order (see Fig. [3]). The 
pairing correlation function -F'a.b(R') rapidly vanishes be- 



yond the nearest neighbor. We also determined the low 
energy contribution of the Drude weight [l^ and found 
-Dorude — S.lSxt for Small X. Despite the addition of mo- 
bile charge carriers, singlets form on diagonal bonds and 
the system remains in a non-superconducting but metal- 
lic phase with strong valence bond correlations. The ab- 
sence of superconductivity for t' > is consistent with 
the finding of Ref.t20;]. 

Correlated Superconductor: The hole doped case, on 
the other hand, shows superconductivity (see Fig. 3). The 
pairing gap has dx2-y2 symmetry with A^^^+x = — A^^^+y. 
Also its magnitude is larger for hole doping compared to 
the electron doped case (we find A^^^+x — t for hole dop- 
ing X = 0.1). The magnetic correlations in this super- 
conductor are very different from the valence bond crys- 
tal. As shown in the inset of Fig.jSj we find (Si • Sj) < 
with rather comparable magnitude for all neighbors, even 
when they are coupled via the exchange interaction J. 

As pointed out above, the origin for the strong asym- 
metry between hole and electron doping is due to the 
constructive {f > 0) and destructive {f < 0) interfer- 
ence of hopping paths in the system. This leads to a 
delocalization of singlets and a true RVB state for hole 
doping while for electron doping the underlying valence 
bond crystal is only moderately affected. 

Inhomogeneous Superconductor: So far we assumed 
that the system is fully homogeneous, i.e. that no spatial 
symmetry in the problem is broken. Next we allow for 
the local potentials xu i effective hopping elements Xij as 
well as pairing gaps A^ to vary in space. Specifically, 
we allow for arbitrary values of these quantities inside 
the crystallographic unit cell, consisting of four Cu-sites. 
For electron doping no additional symmetry is broken. In 
case of hole doping, we find that the spontaneous emer- 
gence of inhomogeneous solutions leads to a reduction of 
the ground state energy {AE ~ 4.3 x 10~^t for x = 0.1 
which is three times the energy gain of the homogeneous 
d-wave state compared to the Ay = state). While the 
charge is homogeneous, we find large variations of the 
pairing gaps following the pattern indicated in the right 
panel of Fig[TJ For x ~ 0.1 we find for strong bonds a gap 
amplitude A ~ 2.3t compared to ~ 0.2t for weak bonds. 
The plaquette state with periodic variation of the pairing 
strength is driven by a gain in the kinetic energy. 

In this plaquette state, the effective hybridizations x' 
on the diagonals also varies in space, x' on diagonals 
surrounded by large pairing gaps is reduced compared 
to x' on plaquettes with small gaps. At half filling, an 
analysis of the Shastry-Sutherland lattice with Sp{N) 
symmetry[13j] led to several plaquette states, including 
one similar to ours, once strong quantum fluctuations 
were taken into account. Doping with charge carriers is 
clearly one route to enhance quantum fluctuations. Phys- 
ically, d-wave RVB fluctuations between the states |=) 
and |ii) inside a plaquette seem strongest if the system 
is allowed to break up in a spatially inhomogeneous pat- 
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tern. 

The most interesting aspect of this new plaquette or- 
dered state is that long range superconducting order ips 
is strongly enhanced in the plaquette state compared to 
the homogeneous state as seen in Fig. [3] The pairing cor- 
relation function i^a.b (R-) reaches its long distant limit 
after several lattice constants. Strong pairing in plaque- 
ttes separated by less that this superconducting coher- 
ence length is therefore very efficient to boost a supercon- 
ducting state For hole doping the magnitude \'ipsf is 
approximately 10% of the value one obtains within BCS- 
theory for the same pairing gap. For electron doping 
the effect of quantum fluctuations much more dramatic 
since \ips\^ = even though a BCS theory with finite gap 
would always yield a finite order parameter. 

A simple analytic approach that provides useful insight 
in case of other doped Mott insulators is the slave boson 
mean field theory. The approach is closely related to 
the wave function, Eqlll with the key difference that the 
projection P is done only on the average. Performing a 
slave boson calculation, we find, in agreement with ear- 
lier work 1J| , that the strong asymmetry between elec- 
tron and hole doping does not emerge. The plaquette 
state found in this paper does not follow from the slave 
boson calculation either. Finally, by using the full projec- 
tion we searched unsuccessfully for a staggered flux phase 
with Xij — Xe^*^ , found in a slave bosons calculation 22 1. 
These results demonstrate that proper projection, which 
enters the slave boson theory via strong gauge field fluc- 
tuations, is crucial for the Shastry-Sutherland lattice. 

In summary, we have identified a doped Mott insulator 
with a very rich behavior as function of varying charge 
carrier concentration. Electron doped systems display 
metallic phases with strong valence bond correlations and 
are the first example for a doped Mott insulator where an 
RVB-wave function yields local pairing but no supercon- 
ductivity. Hole doped systems enter a d-wave RVB state, 
become superconducting and an inhomogeneous bond or- 
der emerges spontaneously which enhances superconduc- 
tivity. From the energy gain of the paired state and the 
values of the exchange energies Jij we estimate a super- 
conducting transition temperature Tc of several Kelvin. 
Experimentally, attempts to dope SrCu2(B03)2 have not 
led to metallic behavior 23]. However, recent advances 
in the sample preparation of transition metal oxides 24 
give us every reason to be optimistic that the obstacles 
to dope this material can be overcome. We hope that our 
results stimulate research in this direction. 
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